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We study surface nanobubbles using molecular dynamics simulation of ternary (gas, liquid, solid) 
systems of Lennard- Jones fluids. They form for sufficiently low gas solubility in the liquid, i.e., for 
large relative gas concentration. For strong enough gas-solid attraction, the surface nanobubble is 
sitting on a gas layer, which forms in between the liquid and the solid. This gas layer is the reason for 
the universality of the contact angle, which we calculate from the microscopic parameters. Under the 
present equilibrium conditions the nanobubbles dissolve within less of a microsecond, consistent with 
the view that the experimentally found nanobubbles are stabilized by a nonequilibrium mechanism. 



When liquid conies into contact with a solid, 
nanoscopic gaseous bubbles can form at the interface: 
surface nanobubbles [JQ- These bubbles were dis- 
covered about 15 years ago, after Parker et al. pre- 
dicted their existence to explain the long-ranged attrac- 
tion between hydrophobic surfaces in water [4j. Many 
Atomic Force Microscopy (AFM) and spectroscopy mea- 
surements have since then confirmed the existence of 
spherical cap-shaped, gaseous bubbles at the liquid-solid 
interface. 

Various open questions remain about surface nanobub- 
bles, and in this Letter we will address three crucial ones: 

(i) How do surface nanobubbles form? This question 
is difficult to answer by experimental means, since the 
formation process is too fast to be observed by AFM. 

(ii) A second question regards the contact angle of sur- 
face nanobubbles which is found to disagree with Young's 
law: all recorded nanobubbles have a much lower gas-side 
contact angle than expected, and seem to be universal 
within 20 degrees, (iii) Finally, AFM showed that sur- 
face nanobubbles can be stable for hours or even days, 
whereas the pressure inside these bubbles due to their 
small radius of curvature (R c ~ 100 nm) would be several 
atmospheres due to the Laplace pressure: Ap = 2^/R c , 
with 7 the liquid-vapour surface tension. A simple cal- 
culation then shows that surface nanobubbles should dis- 
solve within microseconds, which is 9 to 10 orders of mag- 
nitudes off with respect to the experimental data. 

In this paper, we use Molecular Dynamics simulations 
(MD) to study surface nanobubbles in simple fluids. Us- 
ing MD, we are able to answer questions (i) and (ii), and 
provide important information with respect to question 

(iii) . MD arc well-suited for nanobubbles, because the 
temporal resolution is of order fs, and since all atom's 
motions are resolved, the spatial resolution is intrinsi- 
cally high enough to resolve nanobubbles. This atomistic 
model allows to study microscopic details that are inac- 
cessible by experimental means and standard continuum 
mechanics. FigureQ]shows how surface nanobubbles form 
in a typical simulation of a liquid containing gas. The 
gas will homogeneously nucleate to form a bubble, which 



subsequently attaches to the wall. We will analyze the 
nucleation process in detail and quantify how the con- 
tact angle of the bubble changes upon varying gas solu- 
bility. The enhanced gas concentration ( "gas-enrichment 
layer" ) at the solid-liquid interface, which is strongest at 
hydrophobic substrates, will turn out to play a key role 
to account for the universality of the contact angle. 

Numerical setup. The studies in this paper are per- 
formed using simple fluids, which contain no molecules 
but rather separate atoms that interact with each other 
through the Lennard- Jones (LJ12-6) potential: 

"-M(?)"-(?)1- (1) 

Here, e,j and cry are the interaction strength and range 
between particles i and j, respectively. All simulations 
were performed using the GROMACS-software package 
and were done at constant temperature, volume and 
number of particles (T, V, and N constant). The aug- 
mented Berendsen thermostat described in ref. [B| was 
used in all simulations. We verified that this thermo- 
stat yields the same result as the Nose-Hoover thermo- 
stat [6[. In all simulations three types of particles were 
used: fluid I, fluid II, and solid particles. The fluid par- 
ticles (I & II) behave like ordinary particles in a MD- 
simulation and thus move around the system. Contrarily, 
the solid particles are constrained to their initial position 
throughout time and constitute the immobile substrate. 
The interaction parameters of the fluids are chosen such 
that at the temperature considered (T = 300 K) fluid 
I is in the liquid state and fluid II in the gas state, 
and they will be referred to by these states through- 
out the rest of the paper. These interaction parameters 
are: (a ss ,a u ,a gg ) = (0.34,0.34,0.50) nm, (e ss ,e u ,e gg ) = 
(1.2, 1.2, 0.4) AsT, with ks Boltzmann's constant. For 
cross-interactions we use: Uij = (an + (Xjj)/2, and 
(e s i,e sg ,ei g ) = (0.8, 0.8, 0.7)fcsT, unless otherwise stated 
in the text. The cut-off length of the potential function 
was set at r c = 5cr = 1.7nm. 

The time step for the simulations is dt = t/400 « 2 fs, 
where t is a characteristic timescale of atomic motion 
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FIG. 1. (Color) a) Initial conditions: a liquid layer (blue) 
is placed on top of a solid substrate (bottom, red). Gas 
is dissolved inside the liquid layer (green), b) Simulation 
{eig/eu = 0.58, ai g /an = 1.32) after about 0.1 ns: nucle- 
ation occurs, c) t = 10 ns: a surface nanobubble has formed, 
d) Parameter space where the solubility of the gas was tuned 
through the parameters ti g /tu and oigjou. As the gas be- 
comes increasingly soluble (going up, left in Fig. []J) a sudden 
transition takes place where no nanobubbles nucleate. The 
gas then remains in a dissolved state and partially escapes 
to the gas-phase above the liquid layer until equilibrium is 
reached. 



t = cj/ym/eH, with m the mass of the liquid particles 
(20 a.m.u.). The initial conditions are shown in Fig. [Ip,: 
on top of the substrate we place a layer of liquid with dis- 
solved gas. Periodic boundary conditions are present in 
all directions (a;, y, z); the resulting nanobubbles are ap- 
proximately 20-30 nm wide, and 10-20 nm high, depend- 
ing on the contact angle. The formation and behaviour 



of bubbles was found to be independent of simulation 
box-size, which was set at 40x40x5.5 nm 3 . 

Bubble nucleation. What determines whether 

nanobubbles form? We address this question by vary- 
ing the relative interaction strength and the relative in- 
teraction size. We then explore the parameter space to 
see under what conditions nanobubbles form. The result 
is shown in Fig. [T}i. Decreasing ei g /tu (going down in 
Fig. [T]i) results in a lower solubility of gas in the liquid, 
and since the absolute concentration of gas is kept con- 
stant this means that the liquid becomes more and more 
supersaturated. Eventually, homogeneous nucleation oc- 
curs and a nanobubble forms in the bulk liquid phase, 
which finally attaches to the surface. Increasing ai g /au 
(going right in Fig. [TJi) leads to the same effect: due 
to the increased size of the gas atoms it becomes ener- 
getically less favourable to remain dissolved in the liq- 
uid phase. Eventually, when the gas molecules are large 
enough nanobubbles form due to the supersaturation of 
gas in the liquid. From this, we can conclude that a 
local supersaturation of gas inside the liquid is a possi- 
ble mechanism to generate surface nanobubbles. These 
results are consistent with the experimental findings in 
ref. 0, where it was reported that nanobubble formation 
strongly depends on the (relative) gas concentration in 
the liquid. Although the concentration required to spon- 
taneously form nanobubbles is far greater than the satu- 
ration concentration, we point out that during deposition 
of liquid on a substrate gas can be trapped leading to very 
high local transient concentrations which would not be 
reflected in measurements of the global gas concentration 
in the liquid. In fact, numerous experimental papers have 
pointed out that the method of deposition is of great im- 
portance for achieving surface nanobubbles . We have 
to note, however, that other mechanisms not considered 
here can also induce the formation of nanobubbles (e.g. 
heterogeneous nucleation, bulk desorption of gas from 
micropancakes 0, ) ■ The nanobubbles produced in our 
simulations are found to be reproducible and, at the very 
least, can be studied regarding their shape and stability. 

Universal contact angle. Now that we can simulate 
nanobubbles, we focus on the universal contact angle of 
surface nanobubbles found in experiments [UO]. For 
this we use similar initial conditions as in Fig. [TJi, with 
a region in the liquid with a very high gas concentra- 
tion, providing us with control over where the nanobub- 
bles form, and how much gas they contain (about F000 
atoms). We measure the gas-side contact angle[14j at 
varying values of the solid-gas interactions e sg . As can 
be seen in Fig. [5] (triangles), we observe that the gas- 
side contact angle of the nanobubble decreases (i.e. the 
nanobubble becomes flatter) , as the solid-gas interaction 
is increased. This trend saturates around e sg /e/; w 1, 
where the gas-solid attraction matched the liquid-liquid 
attraction [16j. The observed saturation contact angle 
(75°) is close to the contact angle of nanobubbles that is 
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FIG. 2. (Color) The effect of an enhanced gas-solid interaction 
strength. As the interaction strength e sg /eu is increased, the 
adsorbed gas density (pg l , blue squares) increases as well until 
a saturation limit. As the adsorbed gas density increases, the 
gas-side contact angle 6 lowers (cos 6 indicated by red trian- 
gles). The red solid line is a fit to the mean-field expression 
{2} taking into account the screening of the adsorbed gas (see 
text). 




Collapse 
"6 time (ns) 200 



20 40 60 

time (ns) 



80 



FIG. 3. Number of gas molecules n gas inside the nanobubble 
as function of time for nanobubbles on different substrates 
with different equilibrium contact angles. Initially, the fluid 
is supersaturated and a bubble quickly forms within a few ns. 
Shortly after the bubble has formed, the liquid is still su- 
persaturated causing gas to enter the bubble ("Start-up"). 
After about 20 ns the gas in the liquid achieves the equi- 
librium concentration, and the nanobubbles start to dissolve 
("Dissolution"). The dissolution rate of the nanobubbles is 
independent of the contact angle. The inset shows the full 
dissolution of the 6 eq = 93° bubble after 0.2 fis. 



found in experiments (60°). On the same figure we show 
the evolution of the 2D number density of gas concentra- 
tion at the wall inside the liquid phase, p s g . Remarkably, 
this concentration exhibits a trend that is very similar 
to that of the contact angle. Stronger solid-gas interac- 
tions lead to a planar area of high gas concentration at 
the solid-liquid interface, which is called a gas-enrichment 
layer and which has been observed to exist in both sim- 
ple liquids as well as real liquids (l7Hl"9j . In experiments 
high-density gas adsorbates (micropancakes) have also 
been observed 0, Q • Figure [5] shows that the adsorbate 
density increases with e sg , until it finally saturates to a 
2D number density of p s g l = 0.7 atoms per nm 2 in the 
first gas layer above the substrate. 

The increase of the gas density near the wall is indeed 
the origin of the flattening of the nanobubbles. Namely, 
the presence of the gas weakens the attractive interaction 
between solid and liquid molecules: the liquid does not 
'feel' a solid half-space anymore, but there is now a dense 
gas-layer that effectively renders the wall more hydropho- 
bic. This effect can be quantified using the approximate 
equation for the contact angle 



cos 6 g = 1 — 2 



Ps^sl 
PlCll 



(2) 



which can be obtained from a mean field argument. This 



expression contains only the solid and liquid densities p s , 
pi, and the solid-liquid and liquid- liquid interactions e s ;, 
tu. The vapour phase has a negligible contribution (p v 
is small compared to pi and p s ) and the solid-solid in- 
teraction is irrelevant since the solid is non-dcformable. 
In the case of a dense adsorbate, the attraction e s i is re- 
duced to an effective interaction with e s i < e s /. In 
addition, the adsorbate density is lower than the origi- 
nal solid density, due to the large size of the gas atoms, 
and gives an 'effective' density p s < p s . According to 
([2|), both effects lead to a lower gas-side contact angle. 
The solid line in Fig. [2] shows the predicted contact an- 
gle by this expression, assuming an average interaction 
strength: e s i = (e s i + e-ig)/2. Here the effective density is 

estimated by p s = (1 — £r)p s +Pg, as a phcnomenological 
description for the screening of the solid as the adsorbate 
layer density grows. Note that the influence of the vapour 
phase is neglected, as was the case in the model without 
the presence of a gas adsorbate. The model quantita- 
tively explains the observations in MD, in particular, the 
saturation of the contact angle occurs exactly when a 
complete layer of gas adsorbate is formed. It therefore 
provides a very natural explanation for the observed uni- 
versal contact angles in experiments. HUES! 
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Stability. Another aspect of nanobubbles that can 
now be studied is their stability. Are Lennard- Jones 
nanobubbles stable? After formation of the nanobub- 
ble, we use a shape tracker to follow the dynamics of the 
nanobubble. The shape tracker locates a nanobubble by 
performing a circular fit through the liquid-vapour in- 
terface of the curved bubble wall. Different quantities 
can then be computed, such as the radius of curvature, 
the contact angle, the amount of gas inside a nanobub- 
ble, and the angular dependence of gas flux through the 
bubble wall. A good indicator for nanobubble stability 
is the gas content inside the bubble: when the amount 
of gas remains constant the bubble is considered stable. 
The gas contents of nanobubbles on different substrates 
as a function of time are plotted in Fig. |3k. We see that 
none of the nanobubbles are stable; they dissolve on a 
timescale (/is) much shorter than that observed in exper- 
iments (days), see also the inset in Fig. [3] However, this 
fast decay is in agreement with simple macroscopic diffu- 
sion calculations. Furthermore, we find that the contact 
angle of the nanobubbles does not change significantly 
throughout the dissolution process. 

Although the nanobubbles are not stable, some inter- 
esting events occur near the contact line. When studying 
the time-averaged local flux as a function of angular po- 
sition 4> (Fig. HJ}) we see that the flux is highest near 
the contact line, indicating that the substrate plays an 
important role in nanobubble stability. This strong lo- 
calized flux near the contact line is heavily influenced by 
the presence and strength of the gas-enrichment layer, 
which is a plane at the solid-liquid interface in which gas 
atoms can move relatively easily due to a liquid deple- 
tion layer that exists at the same position. The influx 
indicates that there may exist a condition where a dy- 
namic equilibrium is achieved, i.e. the diffusive outflux 
is balanced by the influx at the contact line, explaining 
why in nonequilibrium surface nanobubbles can be dy- 
namically stable [24| (and bulk nanobubbles cannot). A 
coarse exploration of the parameters e s i and e sg has been 
performed in this study, but under the present equilib- 
rium conditions a stable nanobubble was not achieved. 

Of course, there are many more parameters that need 
to be explored, such as the initial radius of the bubble: 
it is possible that nanobubbles below a certain critical 
size are unstable. Also, Lennard- Jones fluids might not 
contain the necessary properties to form stable nanobub- 
bles, such as electrostatic effects. Most importantly, for 
the dynamic equilibrium theory to be true, some driv- 
ing force must exist to sustain the circulation of gas. 
This means that the equilibrium simulations in this study 
need to be adapted to contain such a driving force. Such 
non-equilibrium effects include the presence of a thermal 
gradient (which are likely to be present in experimental 
setups as well) or the formation of gas at the substrate 
(which has been studied using electrolysis [25 - 27 1 ) . 




FIG. 4. Local flux of gas through the liquid-vapour interface 
of a nanobubble attached to a substrate. The gray line indi- 
cates the time-averaged local gas flow direction and strength. 
The gas flow is directed outwards everywhere (iii), except for 
a small region near the contact line [(i) and (ii)] where in a 
very small region a strong in- and outflux are observed, indi- 
cating that there exists a recirculation current. The net effect 
of this recirculation current is found to be of the same order 
as the diffusive outflux. 



Outlook. In conclusion, we have generated and anal- 
ysed the formation and stability of surface nanobub- 
bles in simple fluids. We found that in heavily gas- 
supersaturated liquids nanobubbles nucleate sponta- 
neously which can then migrate towards the surface. In 
experiments, when water is deposited on the substrate, 
it is possible that some gas becomes trapped near the 
solid-liquid interface leading to the required supersatu- 
ration. Other formation mechanisms that cannot be ac- 
cessed by MD can however not be excluded, further work 
is required on this question. The universal contact an- 
gles that surface nanobubbles exhibit in experiments can 
be explained by a dense layer of gas at the solid-liquid 
interface, which has been shown to exist for real liquids, 
that effectively alters the substrate chemistry. Although 
the Lennard- Jones nanobubbles are unstable, some inter- 
esting local gas flows are present near the contact line. 
These gas flows are caused by the symmetry breaking 
due to the solid substrate, and hint towards a dynamic 
equilibrium condition where the diffusive outflux is com- 
pensated by an influx near the contact line. Since an 
energy input is required to sustain a circulatory gas flow 
suggested in the dynamic equilibrium theory by Brenner 
and Lohse [23 |. it is likely that stable nanobubbles can 
only occur in non-equilibrium systems. Simulations of 
non-equilibrium systems, and of systems containing re- 
alistic fluids must be performed to address the question 
regarding the long lifetime of surface nanobubbles. 
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